Spin diode behavior in transport through single-molecule magnets 
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Abstract. - We study transport properties of a single- molecule magnet (SMM) weakly coupled 
to one nonmagnetic and one ferromagnetic lead. Using the diagrammatic technique in real time, 
we calculate transport in the sequential and cotunneling regimes for both ferromagnetic and 
antiferromagnetic exchange coupling between the molecule's LUMO level and the core spin. We 
show that the current flowing through the system is asymmetric with respect to the bias reversal, 
being strongly suppressed for particular bias polarizations. Thus, the considered system presents 
a prototype of a SMM spin diode. In addition, we also show that the functionality of such a device 
can be tuned by changing the position of the molecule's LUMO level and strongly depends on the 
type of exchange interaction. 



Introduction. — Due to their particular physical 
properties, such as an energy barrier for the spin re- 
versal or long relaxation times [1 , single- molecule mag- 
' nets (SMMs) are inherently predestined for applications in 
novel molecular electronic or spintronic circuits [2j[3]. Up 
to now, several different physical mechanisms employing 
SMMs as a key component have been theoretically consid- 
ered. It has been in particular shown that the SMM's spin 
- can be reversed by means of spin polarized currents [4-6 , 
' or by applying a spin bias [7 . Furthermore, when bridged 
between two nonmagnetic metallic leads a SMM can work 
_ as a spin filter [8"^ "9^. Recent experiments on electronic 
< transport through SMMs connected to metallic but non- 
magnetic leads [TOHIl] clearly show that the aforemen- 
tioned ideas are in principle experimentally feasible. From 
a conceptual point of view, the simplest realization of a 
system consisting of a SMM attached to one ferromagnetic 
and one nonmagnetic reservoir could be a device involving 
the scanning tunneling microscope (STM) with a magnetic 
tip and a SMM on a metallic but nonmagnetic substrate, 
Fig. [U^a). The advantage of such a geometry is that by 
choosing an appropriate ligand shell for the molecule, one 
can obtain the specific orientation (e.g. parallel) of the 
molecule's easy axis with respect to the surface [TF. Fur- 
thermore, there are experimental techniques which allow 



for deposition of a film of well-dispersed SMMs on a sub- 
strate, enabling access to individual molecules with the 
STM tip [ 16H2Qj . Thus, it is interesting to consider how 
transport properties of a SMM change if one of two metal- 
lic nonmagnetic leads is replaced by a magnetic one. 

In this Letter we study in general transport through 
a SMM weakly coupled to electrodes with unequal spin 
polarizations. As already shown in the case of quantum 
dots pTH2^ , transport properties of such systems exhibit 
a significant asymmetry with respect to the bias reversal. 
In addition, due to coupling to ferromagnetic leads and 
the spin dependence of tunneling processes, the current 
flowing through such a diode becomes spin polarized and, 
interestingly, the spin polarization may change with re- 
versing the bias voltage. In fact, very recently spin diode 
behavior was predicted and observed experimentally in 
another class of molecular structures, namely in single- 
wall carbon nanotubes [25l|26]. In this work we propose 
a molecular spin diode based on single molecule magnets. 
The considered system consists in particular of a SMM 
weakly coupled to one nonmagnetic and one ferromagnetic 
lead of high spin polarization; Fig.fljb). To determine the 
transport properties of such a system we employ the real- 
time diagrammatic technique (RTDT), which allows us to 
systematically take into account the sequential and cotun- 
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The molecule is thus fully characterized by Hamiltonian 



Metallic nonmagnetic substrate 




Fig. 1: (Color online) (a) Illustration showing a possible real- 
ization of a SMM-based device employing one metallic mag- 
netic lead: the SMM molecule attached to a metallic nonmag- 
netic surface is pinned from the top by a magnetic STM tip. 
For simplicity, the magnetization of the tip and the molecule's 
easy axis are assumed here to be parallel, (b) Schematic of 
a SMM attached to one nonmagnetic (left) and one ferromag- 
netic (right) lead. The current flows due to tunneling through 
the LUMO level which is exchange coupled to the molecule's 
core spin S. 



neling processes contributing to current. The sequential 
tunneling dominates the current for voltages larger than 
the Coulomb correlation energy and becomes exponen- 
tially suppressed in the Coulomb blockade regime, where 
transport is mainly due to cotunneling processes. As pre- 
sented in the following, unequal couplings to the leads give 
rise to a pronounced asymmetry of spin-polarized current 
with respect to the bias reversal. Furthermore, in addition 
to the analysis of system's I-V transport characteristics, 
we also discuss the zero-frequency noise of tunneling cur- 
rent associated with discreteness of charge carriers (shot 
noise). The shot noise proves to be a source of useful 
information about various transport properties of the sys- 
tem, such as effective charges, coupling strengths or vari- 
ous types of correlations, which is hardly attainable from 
direct measurement of electric current ^ . 



Model and theoretical method. — We assume 
that transport through the molecule occurs only via the 
lowest unoccupied molecular orbital (LUMO) level, de- 
scribed by the energy e and the local spin operator s = 
\ Xlcrcr' ^"^(T^crcr'Ccr' 1 whlch Is couplcd to the molcculc's in- 
ternal spin S via exchange interaction J. In principle, the 
exchange interaction can be either ferromagnetic or anti- 
ferromagnetic, depending on the specific type of a SMM. 
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where D is the uniaxial anisotropy constant of a neu- 
tral molecule, Di and D2 stand for corrections to the 
anisotropy due to single and double occupation of the 
LUMO level, respectively, while U describes the Coulomb 
correlations between two electrons occupying the LUMO 
level. We note that although the corrections Di and 
D2 to the anisotropy constant D are not crucial in ob- 
serving the presented effects, we take them into account 
to make the model realistic. In Eq. ([!]) we neglected 
small transverse anisotropy (61. Furthermore, it is as- 
sumed that the easy axis of the SMM is collinear with 
the magnetization of the right electrode, see Fig. [Hb), 
and both electrodes are described by noninteracting elec- 
trons, Hei = EgEk,a4a ^l^a' ^ith being the 
energy of an electron with the wave vector k and spin 
<j belonging to the lead q. Finally, the Hamiltonian de- 
scribing tunneling processes between the molecule and 
the leads reads, = Ek,. + ^g4<], 

where Tq denotes the tunnel matrix elements between the 
molecule and the gth lead. Coupling between the molecule 
and each lead can be expressed by F^ = 27r|Tgpp^ and 
Lg = (F^ + Fi)/2. The spin polarization of the right mag- 
netic lead is defined as p = {p^ — p^)/{p^ + P-)^ where 
is the density of states for the majority (minority) 
electrons. 

To calculate the current flowing through the system 
we employ the real-time diagrammatic technique [271128]. 
This approach allows us to take into account the first and 
second order tunneling processes through SMM in a fully 
systematic way [29]. Within the RTDT, the occupation 
probabilities of finding the system in a many-body 
state Ix) are given by. 



(2) 



where X) is the self-energy matrix accounting for various 
tunneling processes in the system, with a row xo replaced 
by (F, . . . , F) due to normalization of probabilities. The 
current can be then calculated as 



(3) 



with I]^ being the modified self-energy to take into account 
the number of electrons transferred through the system. 
Also the expression for the current noise in the limit of 
low frequencies [30] . 
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dt[{mm^i{o)m)-2{i)' 



(4) 



can be conveniently reformulated in the language of 
RTDT, for details see Ref. [31 . At this point, we ought to 
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Fig. 3: (Color online) The zth component of the SMM's total 
spin St — Sz-\- (c|c^ — c|c^)/2 as a function of the bias voltage 
V and the LUMO level position e for (a) ferromagnetic and (b) 
antiferromagnetic exchange coupling. Hollow arrows depict the 
preferred direction of the SMM's spin with respect to the right 
lead's magnetic moment, whereas the dashed lines correspond 
to the position of main peaks in the differential conductance, 
Fig. [2Ub,d) - i.e. they divide regions representing different 
occupation Q of the LUMO level. The remaining parameters 
are as in Fig. [JJ 



Fig. 2: (Color online) (a,c) The absolute value of the current 
/ in units of Iq = 2eV /h ^ 0.5 nA, and (b,d) differential con- 
ductance G as a function of the bias voltage V and the LUMO 
level position e for (a)-(b) ferromagnetic (J > 0) and (c)-(d) 
antiferromagnetic (J < 0) exchange coupling, \J\ = 0.2 meV. 
Q in (a,c) represents the average charge accumulated in the 
LUMO level. The other parameters are: = 2, D = 50 /ieV, 
Di -5 /ieV, D2 = 2 /xeV, U 1 meV, /cbT 40 /xeV, 
p = 0.9, and r = Fl = Fh = 1 /xeV. 

note that in the next section, rather than the shot noise 
we present the deviation of the current noise from its Pois- 
sonian value, described by the Fano factor F = S/2e\I\. 

Finally, to determine the current order by order in tun- 
neling processes, one performs perturbation expansion of 
the self-energies and probabilities in the coupling strength 
F [271EH]- Here, we have taken into account the first and 
second order terms of expansion, which correspond to se- 
quential and cotunneling processes, respectively [29,32 . 

Results and discussion. — Since the model un- 
der consideration applies to various types of SMMs, here 
we present results only for the conceptually easiest case, 
though capturing essential physics, i.e. for a molecule with 
S = 2 and strong uniaxial magnetic anisotropy. Moreover, 
to keep the discussion most intuitive, we neglect the role 
of electron charge sign, assuming e > 0. 

First of all, we note that due to large spin asymmetry in 
coupling of the SMM to the ferromagnetic lead, Fig.[T]^b), 
the tunneling probability for spin-majority (spin- up) elec- 
trons is much larger than that for spin-minority (spin- 
down) electrons. On the other hand, the rate of tunneling 
processes between the molecule and the nonmagnetic lead 



is the same for both spin orientations. This generally leads 
to an asymmetry of tunneling current with respect to the 
bias reversal. 

In Fig. [2fa)-(b) we display the current and differential 
conductance as a function of the LUMO level position and 
bias voltage for ferromagnetic ( J > 0) interaction between 
the LUMO level and the SMM's spin. Close to the first res- 
onance, e ~ £01 7 see Fig.[2fb), current can flow easily from 
the ferromagnetic lead to the nonmagnetic one {V < 0), 
while it is suppressed for the opposite direction (V > 0); 
see also solid line in Fig. HJ^a). To understand this behav- 
ior one should take into account the following facts. First, 
the molecule's states with one extra electron in the LUMO 
level correspond to the total spin number S = 5/2 and 
S = 3/2, with the former being of lower energy. Second, 
orientation of the molecule's spin depends significantly on 
the current direction, Fig.[3fa), and for V > the SMM's 
spin tends towards antiparallel orientation with respect 
to the electrode's magnetic moment, whereas for V < 
the spin prefers the parallel alignment. Thus, for low bias 
transport can occur mainly via the state corresponding to 
S = 5/2, and this requires spin-down electrons for V > 
and spin- up electrons for V < 0. Unfortunately, tunnel- 
ing rate for spin-down electrons is significantly reduced 
due to fewer available states in the minority spin band 
of the ferromagnetic lead, which effectively leads to the 
suppression of current for ^ > 0. When bias voltage in- 
creases, then the state corresponding to 5* = 3/2 becomes 
active in transport as well, leading only to a small increase 
of the current due to the spin selection rules for tunnel- 
ing processes. However, the blockade for positive bias is 
removed when double occupancy of the LUMO level is ad- 
mitted, which takes place for bias voltages exceeding some 
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Fig. 4: (Color online) Selected cross-sections of density plots in Fig. [2] for specific values of the LUMO level position s, which 
depict the current / flowing through the system (a,d), differential conductance G (b,e), and Fano factor F (c,f) as a function 
of the bias voltage V for the ferromagnetic (left panel) and antiferromagnetic (right panel) exchange couplings. The other 
parameters are the same as in Fig. [2] At low bias voltages the shot noise is dominated by thermal fluctuations while the current 
vanishes leading to a divergency in the Fano factor when V ^ 0. 



threshold value. Furthermore, the electron flow from fer- 
romagnetic lead to the nonmagnetic one is spin polarized 
and degree of this polarization depends mainly on the spin 
polarization of the lead. In fact, in the case of a perfect 
halfmetallic ferromagnet {p — > 1), the current flowing to- 
wards halfmetallic lead would be totally blocked. More- 
over, different time scales associated with spin-majority 
and spin-minority electrons lead in turn to considerable 
current fluctuations and super-Poissonian shot noise, see 
Fig. ijc) for e = 0.2 meV. 

The situation becomes significantly different when the 
LUMO level is doubly occupied in equilibrium; e = —1.5 
meV in Fig. IH^a). Now, the behavior of current is re- 
versed as compared to the case of Q = 0, since the current 
is suppressed for F < 0, i.e. for electrons tunneling from 
the magnetic lead. This is associated with the fact that 
an electron first has to tunnel out of the LUMO level and 
then another electron can enter the molecule. Thus, for 
positive bias a spin-up electron can easily tunnel out to 
the ferromagnetic lead. On the other hand, when the bias 
is reversed and the spin-down electron tunnels out of the 
molecule leaving it in the state corresponding to S = 5/2, 
the current becomes suppressed, as the rate for tunnel- 
ing of spin-down electrons from the ferromagnetic lead to 
the molecule is relatively small. This also leads to super- 
Poissonian shot noise, as shown in Fig. Hl^c). 

More complex transport characteristics are observed 



when the exchange interaction is antiferromagnetic (J < 
0); Fig. [2fc)-(d). The most striking difference is the 
appearance of additional peaks in the current when the 
LUMO level is initially either empty or doubly occupied. 
Figs. [2fc) andlD^d), which are accompanied by negative 
differential conductance (NDC); Figs.EJd) andlU^e). Con- 
sider first the case of empty LUMO level in equilibrium, 
e = 0.2 meV in Fig. lU^d). The key difference is that 
now the molecule's state corresponding to the total spin 
number S = 3/2 has lower energy and determines trans- 
port properties at low voltages. Thus spin-up electrons 
are involved in charge transport for V > and spin- 
down electrons for V < 0. Consequently, the current is 
suppressed for negative bias and can easily flow for posi- 
tive one. When the bias voltage reaches values admitting 
transport through the S = 5/2 state, the current for posi- 
tive bias becomes suppressed by a spin-down electron tun- 
neling to the LUMO level, while suppression for negative 
voltage becomes then lifted. In turn, when bias increases 
further admitting doubly occupation of the LUMO level, 
the blockade for positive bias becomes removed as well. 
Transport characteristics for doubly occupied LUMO level 
in equilibrium can be explained in a similar way. 

On the other hand, when the LUMO level is singly 
occupied and its position corresponds approximately to 
the middle of the Coulomb blockade [dashed line in 
Figs. |4l^a,d)], the diode behavior disappears and the cur- 
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rent recovers symmetry with respect to the bias reversal 
regardless of the type of the exchange interaction J. This 
is due to the fact that now with increasing the bias volt- 
age all charge states of the molecule start taking part in 
transport at the same time, i.e. once the bias voltage 
reaches the threshold. The transport characteristics be- 
come then symmetric with respect to the bias reversal and 
the noise is rather sub-Poissonian, indicating the role of 
single-electron charging effects in transport. Note, how- 
ever, that in the Coulomb blockade regime bunching of 
inelastic cotunneling processes may still result in enhance- 
ment of the shot noise; see dashed lines in Figs.|4l^c,f). Fur- 
thermore, it is visible that in the case of e = —0.5 meV 
the enhancement is much more pronounced in the case of 
the antiferromagnetic coupling. Fig. IH^f). Such a behav- 
ior stems from the fact that for J < both energetically 
lowest lying molecular magnetic states = ±3/2 allow 
the possibility of occupying the LUMO level by an elec- 
tron either with the spin up or down, whereas for J > 
the state = +5/2 {S^ = —5/2) can only accommo- 
date an electron with spin up (down). For this reason, in 
the former case both majority and minority electrons of 
the ferromagnetic lead can participate in transport, thus 
increasing the fluctuations. 

To conclude, we have proposed a molecular spin diode 
based on single molecular magnets coupled to one ferro- 
magnetic and one nonmagnetic leads. We have shown that 
the spin-dependent current flowing through the device be- 
comes suppressed for one bias polarization and is enhanced 
for the opposite one. The suppressed transport is then ac- 
companied with super-Poissonian shot noise. Moreover, 
we have demonstrated that the transport characteristics 
of SMM spin diodes strongly depend both on the number 
of electrons occupying the LUMO level, and on the type 
of exchange interaction between the LUMO level and the 
SMM's core spin. This gives the possibility to tune the 
functionality of a SMM spin diode by shifting the position 
of the LUMO level with a gate voltage. 

Finally, we would like to note that although the pre- 
sented results were calculated for a specific molecule with 
5* = 2, the spin diode behavior can be observed in a 
large class of molecules. The main requirement is to have 
considerable spin asymmetry in the couplings to the left 
and right electrodes and well-defined spin states in the 
molecule, i.e. J^D ^ T, where T is the experimental 
temperature. 
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